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Determination of the Amplitude-Probability Distribution 
of Atmospheric Radio Noise From Statistical Moments 



W. Q. Crichlow, C. J. Roubique, A. D. Spaulding, and W. M. Beery 

(July 23, 1959) 

During the International Geophysical Year, the National Bureau of Standards estab- 
lished a network of atmospheric noise recording stations throughout the world. The ARN-2 
noise recorder at these stations measures three statistical moments of the noise: average 
power, average voltage, and average logarithm of the voltage. An empirically-derived 
graphical method of obtaining an amplitude-probability distribution from these three 
moments, and its development, is presented. Possible errors, and their magnitudes, are 
discussed. 



1. Introduction 

The interference to reception of radio signals 
caused by atmospheric noise depends not only on 
its average level but also on its detailed characteris- 
tics. A complete description of the noise received 
at a particular location would require an exact 
determination of the variation of its instantaneous 
amplitude as a function of time. Because of the 
complexity of the noise structure and the fact that 
no two samples of noise have identical amplitude- 
time functions, it is necessary to resort to simpler 
types of description of a statistical nature. 

Various methods of measurement have been in- 
vestigated by experimenters throughout the world 
for many years, and the relationship between the 
measured values and the interference caused to 
various types of service has been studied. 

In order to determine the most significant charac- 
teristics to be measured, an international subcom- 
mittee was formed within Commission IV of URSI * 
to study the problem. A recommendation [l], 2 based 
on the subcommittee's report [2], was endorsed by 
Commission IV at the Xllth General Assembly of 
the URSI in 1957. In the Annex to this recommenda- 
tion, it is stated that "it would be highly desirable to 
obtain detailed statistical information on the ampli- 
tude-probability distribution of the instantaneous 
envelope voltage, the various time functions, and the 
direction of arrival of the noise at many locations 
throughout the world. However, since continuous 
detailed measurements of this type at all frequencies 
and a large number of stations become prohibitive 
because of the complexity of the equipment involved 
and the large number of personnel necessary to carry 
out the observations, it is recommended that these 
complete studies of the detailed noise characteristics 
be confined to a few selected locations with continuous 
measurements of one or more simple parameters at a 
large number of stations. " 



i International Scientific Radio Union. 

2 Figures in brackets indicate the literature references at the end of this paper. 



During the International Geophysical Year, the 
National Bureau of Standards established a network 
of 16 noise recording stations distributed throughout 
the world [3]. Measurements are made with the 
ARN-2 recorder which was designed at NBS to 
record the average noise power in a bandwidth of 
about 200 cps on eight fixed frequencies between 
13 kc and 20 Mc. In order to obtain additional 
information on the character of the noise, two other 
statistical moments, the average envelope voltage, 
and the average logarithm of the envelope voltage 
are recorded at ten of the stations. These data are 
available from IGY World Data Center A, NBS, 
Boulder, Colo., at a nominal fee. The results of 
data analysis will also be published in an NBS 
Technical Note for public sale. 

The amplitude-probability distribution is very 
useful in assessing the interference potential of noise 
and has been measured at several locations [4, 5, 6, 7]. 
Since data on the moments measured by the ARN-2 
can be recorded continuously over the full frequency 
range with much less equipment and fewer personnel 
than would be required for a comparable amount of 
data on the complete amplitude-probability distribu- 
tion, an investigation was made at the NBS of 
methods of deriving the complete distribution from 
these three moments. 

It is the opinion of the authors that mathematical 
models for the distribution [1,5,7,8] have proved to 
be extremely cumbersome and of doubtful value in 
deriving the distribution from measured moments; 
therefore, an empirically derived graphical method 
has been developed. 

2. Measurements of the Amplitude- 
Probability Distribution 

Simultaneous measurements of the amplitude- 
probability distribution and the three moments 
measured by the ARN-2 were needed over a wide 
frequency range with a variety of types of noise. 
There were two specific purposes in making these 
measurements: 
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a. The distribution measurements, alone, were 
needed to establish typical shapes that could be ex- 
pressed by numerical parameters, easily related to 
the measured statistical moments. 

b. The simultaneous measurements of the distri- 
butions and the moments were needed to check the 
accuracy of the ARN-2 noise recorder in measuring 
the moments. 

Accordingly, a distribution meter was constructed 
to be used in conjunction with the ARN-2 recorder. 

The operation of the distribution meter is based 
on counting techniques. The total number of IF 
cycles in the ARN-2 which exceed a specified ampli- 
tude in a given time interval are counted. By proper 
choice of the time interval, the percentage of time a 
level is exceeded can be read directly from the 
counter. Three counters are used to provide simul- 
taneous measurement at three levels, separated in 
amplitude by 6 db. By use of a calibrated, man- 
ually-operated attenuator, the total possible range 
of 100 db can be covered in levels separated by 
multiples of 2 db. The time duration of each sam- 
ple can be preset to provide an automatically timed 
sample having an interval ranging from 1 sec to 1 hr, 
in 1 sec increments. 

The actual intervals were chosen of sufficient dura- 
tion to provide a statistically significant sample of a 
slowly varying phenomenon such as thunderstorm 
activity. Longer intervals were required for the 
high-level, low-probability noise pulses than for the 
low-level, high-probability noise. The time neces- 
sary to obtain a complete curve ranged from about 
15 min to 30 min. The longer period was necessary 
for measuring noise with a large dynamic range. 

A check for accuracy of operation was made by 
obtaining a distribution of thermal noise. The 
Rayleigh distribution thus obtained confirmed the 
accuracy of measurement. Figure 1 is a typical 
amplitude-probability distribution of atmospheric 
radio noise at 13.3 kc. This was plotted on Rayleigh 
graph paper, whose coordinates were chosen so that 
a Rayleigh distribution will plot as a straight line 
with a slope of — ){. z The coordinates are shown as 
noise level in decibels above the root mean square 
voltage versus the percentage of the time that each 
level is exceeded. This measurement was made 
with the distribution meter at the Boulder Labora- 
tories. The dynamic range between the 0.0001 per- 
cent and 99 percent intercepts is approximately 84 db. 

It was found that these amplitude-probability 
distributions for atmospheric radio noise can be ade- 
quately represented by a three-section curve which 
has been drawn through the data points in figure 1. 
This curve has a characteristic shape that can be 
described by means of four numerical parameters to 
be defined later. The accuracy of fit using this par- 
ticular kind of curve composed of three sections is 
typical of more than 100 distributions experimentally 
measured at Boulder. In the particular case of 
thermal noise, this curve becomes a straight line. 



3 These coordinates are log of voltage versus— }4 logio (—In of probability). 
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PERCENTAGE OF TIME ORDINATE IS EXCEEDED 

Figure 1. Measured amplitude-probability distribution of 
atmospheric radio noise. 

Data obtained by other experimenters [2] in Alaska, 
Panama, England, and Florida when plotted on 
Rayleigh paper also fit curves of this form. 

The lower portion of the curve, representing low 
voltages and high probabilities, is composed of many 
random overlapping events, each containing only a 
small portion of the total energy. Therefore, this 
portion of the curve must approach a Rayleigh 
distribution. 

The section representing very high voltages ex- 
ceeded with low probabilities is, in general, composed 
of nonoverlapping large pulses occurring infrequently. 
From experimental measurements of atmospheric 
noise distributions, this section has been found to be 
well represented by a straight line on Rayleigh graph 
paper. On this graph paper, the remaining section 
of the distribution has been found to correspond 
quite closely to an arc of a circle tangent to the 
above two straight lines. Based on experimental 
evidence, it is assumed that no pulses of amplitude 
greater than those for which the probability^ of being 
exceeded is 10~ 6 are present in the distribution. 

3. Definition of Parameters and Statistical 
Moments 

Since on the coordinate paper used for the plotting 
of amplitude distributions for atmospheric radio noise 
a Rayleigh distribution will always have a slope of 
— ji, and since the center of the circular arc men- 
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tinned above will lie on the bisector of the angle 
Conned by the Rayleigh and high-voltage, low- 
probability lines, four parameters are needed to de- 
line the distribution. A point through which the 
Rayleigh line passes, a point and a slope for the 
high-voltage, low-probability line, and a parameter 
of some kind to determine the radius of the circular 
iwv are needed. 

These four parameters are defined as follows: 

appoint at which the Rayleigh line intersects the 
0.5 probability coordinate. 

b = point on that tangent to the circular arc which 
is perpendicular to the bisector of the angle formed 
by the Rayleigh and high-voltage, low-probability 
lines and oil the probability coordinate through the 
vertex of this angle. 

c=point where the high-voltage, low-probability 
line intersects the 0.01 probability coordinate. 

s = slope of the high-voltage, low-probability line. 

The three statistical moments experimentally 
determined are defined as follows: 



mis voltage: e T 



■=vr 



e 2 dp 



average voltage: e &ve = edp 
antilog of the mean log of voltage: 

= antilog log e dp. 



^log - 



where e is 1/^2 times the instantaneous envelope 
voltage, and p is the probability of e being exceeded. 
Since the ARN-2 noise recorder measures the db 
difference betweeen the e rm3 and e ave and between the 
e rms and e Yog , and in order to minimize the necessary 
graphical calculations, the above parameters are 
altered to relate them to the Rayleigh section of 
the distribution as follows (sec fig. 2): 

^4=20 log — =db difference between Rayleigh at 

^rms 

0.5 probability and the e rms value. 

B=db difference between the point of intersection 
of the Rayleigh and high-voltage, low-proba- 
bility lines and the tangent to the circular arc. 
(The tangent is drawn perpendicular to the 
bisector of the angle formed by the Rayleigh 
and high-voltage, low-probability lines.) 

C=db difference between c and value of Rayleigh 
at the 0.01 probability coordinate. 

X= absolute value of slope of high-voltage, low- 
probabilitv line relative to Rayleigh. 
X=-2s 

The moments actually measured are labeled as 
follows: 

F„=20 log^- 8 



Z<«=20log- 



'log 
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X = 


Absolute volue of slope of 
high voltage, low probability 
line relative to Rayleigh line. 


\ 






























^^ 


\ 
























^ 




/ 




























/ 

/ 




/ 
/ 








I 
















/ 


/ 

/ 
/ 






( 


\ 

\ 








/' 
/ 

A 








/ 

/ 
/ 














Y 






._ 


N 








/ 


/ 

f 














\ 
















/ 
/ 







































-40 
0.0001 0.01 0.1 I 5 10 20 30 40 50 60 70 80 90 95 98 99 

PERCENTAGE OF TIME ORDINATE IS EXCEEDED 

Figure 2. Definition of parameters for amplitude-probability 
distribution atmospheric radio noise. 

The e rma is also measured. 

The following relation will always hold true for the 

distribution: 

therefore, V d and L d are always positive. 

4. Determination of the Parameters as 
Functions of the Measured Moments 

Since four parameters are needed and only three 
statistical moments are measured, a dependency 
between two of the parameters must be determined 
in order to be able to obtain tin 4 remaining three 
independent parameters as functions of the measured 
moments. From experimentally measured distribu- 
tions at eight frequencies and various bandwidths, 
the following linear relation between .V and B is 
found to hold in general (fig. 3) . 

£=1.5(X-1) 

Distributions were employed in which C was 
allowed to vary from to 40 db and X from 1 to 12. 
These distributions were numerically integrated in 
order to determine their three moments. From these 
integrations, V d and L d were obtained as functions 
of X and C. A was also obtained as a function of 
X and O (fig. 4). These functional relations were 
solved simultaneously by graphical means to obtain 
X and C as functions of V d and L d (figs. 5 and 6). 
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Figure 3. Experimental correlation of B and X. 

Figure 7 shows the distribution obtained by the 
above procedure for the following measurements: 

V d = 4.5 db below e Tms 

L d = 9.2 db below e Tms . 

On the X versus V d and L d curves, the above meas- 
urements are seen to give an X of 3.7, and on the C 
versus V d and L d curves the above measurements 
give a C of 15.8 db. Using the A versus X and C 
curves, this pair of values for X and C gives an A of 
12.3 db. Point a is located on the 50 percent coordi- 
nate 12.3 db below the E^ and the Rayleigh line 
with slope — ]i is drawn. Point c is located 16 db 
above the point where the Rayleigh line intersects 
the 1 percent coordinate, and through point c the 
high-voltage, low-probability line with slope s= 
— X/2 = — 1.8 is drawn. The angle formed by the 
Rayleigh line and the high-voltage, low-probability 
line is now bisected and B is calculated from the 
relationship B=l.b(X-\). B is found to be 3.6 db. 
Point b is located 3.6 db above the intersection of the 
Rayleigh and high-voltage, low -probability lines and 
through this point the tangent line is drawn perpen- 
dicular to the angle bisector. Distance 2 to 3 is 
made equal to distance 1 to 2, and at point 3 on the 
Rayleigh line a perpendicular to the Rayleigh line 
is constructed. This perpendicular intersects the 
angle bisector at the center of the circular-arc portion 
of the distribution. The complete distribution is 
now determined with the construction of this arc. 



It should be noted from the curves (figs. 5 and 6) 
that for a given V d only a certain range of values of 
L d will result in a distribution of the above form. 
These ranges are enclosed approximately within the 
dashed lines. Also, the combination of the above 
given V a and the minimum allowable L d results in a 
nonunique solution, i.e., there is, essentially (within 
the range of accuracy of measurements), an infinite 
number of distributions, all with the same Y d and 
L d of the above special combination. Fortunately 
these special combinations almost never occur. 

The above facts allow us to determine to some 
extent the validity of data received from recording 
stations, and whether a measured sample is true 
atmospheric noise or contaminated noise. 

5. Errors in Noise Measurements and Their 
Influence on the Calculated Distribution 

The accuracy with which a distribution can be 
determined from measured moments depends upon 
(1) the validity of the form factor which has been 
taken to represent the distribution, and (2) the 
accuracy of the measured statistical moments. 

The validity of the form factor is very difficult to 
check because of the fact that 15 to 30 min are 
required to measure the distribution, and the sta- 
tistics of the noise do not always remain constant 
for this period of time. As a check for changing 
noise characteristics, each series of distribution 
measurements was immediately repeated and any 
changes noted. Further, the moments recorded by 
the ARN-2 provided information on the stability 
of the noise statistics. Figure 1 represents an 
example of the fit in the case where the noise moments 
were measured precisely and the noise characteristics 
appeared to be constant during the interval required 
for the measurement of the entire distribution. 
Methods of recording noise samples on magnetic 
tape are under development and these samples will 
enable a more critical study of the form factor as 
well as other detailed characteristics. 

Errors in the ARN-2 could be detected by com- 
paring the recorded three parameters with the same 
three parameters as numerically integrated from the 
curve recorded by the distribution meter. The most 
likely source of error in the ARN-2 was the possible 
saturation of the square law detector with a resultant 
lowered value for E Tms . Since the average voltage 
and average logarithm of the voltage are measured 
in decibels relative to E Tms , their values would be in 
error by the same amount. 

Approximately sixty sets of distribution measure- 
ments were analyzed for the error made by the ARN- 
2 recorder in measuring the moments. These 
ARN-2 errors were divided into class-intervals and 
plotted as error-probability curves on arithmetic 
probability graph paper in figure 8. The fairly 
straight lines obtained indicated largely random er- 
rors, with some bias due to square-law detector satu- 
ration. For the E rm8 values, the standard deviation 
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Figure 4. A versus X and C. 



was about 1.1 db, with an average error of —0.4 db. 
Ninety percent of the values were within ±2 db 
of the mean. For the E &ve values, the standard 
deviation was about 0.8 db, with an average error 
of +0.2 db. Ninety percent of the values were 
within ±1.5 db of the mean. For the E log values, 
the standard deviation was about 1.4 db, with an 
average error of +0.6 db. Ninety percent of the 
values were within ±2.5 db of the mean. 

An effort was made to correlate the error with 
dynamic range between the 0.0001 percent and 99 



percent intercepts of the amplitude-probability dis- 
tributions. The same data were divided into class- 
intervals and the medians of the class-intervals which 
contained sufficient samples for statistical signifi- 
cance were plotted against dynamic range in decibels 
in figure 9. The best straight -line fit indicates good 
correlation of error with dynamic range, a negative 
error for E rm8 and a positive error for E &ye and 
E log . The error in E TmB at 85-db dynamic range 
varied from —2.7 to +2.1 with a median value of 
about 1 db. 
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Figures 10 and 11 show the influence of errors in 
the moments measured by the ARN-2 on distribu- 
tions determined from the moments. In these 
figures, the measured moments are indicated as 
fi'rms etc., and the moments obtained from the 
simultaneously measured distribution are indicated 
by the same symbols without the primes. The 
examples chosen contain errors typically measured 
by the ARN-2. In figure 10, the errors in the 
measured moments are small and the distribution 
derived from the slightly erroneous ARN-2 moments 
fits the measured points with negligible error. In 
figure 11, the dynamic range is wider with corre- 
spondingly larger errors in the ARN-2 measured 
moments. Thus, the departure of the distribution 
derived from the ARN-2 measured moments is 
greater. These departures, though significant are 
not considered excessive in view of the fact that 
they are of the same order of magnitude as the 
changes in the noise itself during an hour. 



Thus, it can be concluded that this is a satis- 
factory method of deriving the complete amplitude- 
probability distribution from statistical moments, 
since the departures from the true distribution will 
be small when the moments are measured accurately. 

6. Bandwidth Considerations 

The foregoing method of obtaining an amplitude- 
probability distribution from the three measured 
moments results in a distribution which is only valid 
for the bandwidth in which the moments were 
measured. Fulton [8] developed a method of band- 
width conversion which was accurate for the high- 
amplitude, low-probability portion of the distribu- 
tion, but failed to convert the lower portion of the 
curve as accurately since it was an empirical method 
based upon insufficient measurements. A study is 
now in progress for a more accurate method of 
bandwidth conversion. 
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Figure 7. A?nplitude-probability distribution determined for 
V d =4.5 db and L d = 9.2 db. 



Figure 6. C versus La and V d . 
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Figure 8. Distribution of errors of measured moments. 
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Figure 9. Error of measured moments versus dynamic range. 
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Figure 10. Comparison of measured distribution with that 
obtained from moments measured by ARN-2. 
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Figure 11. Comparison of measured distribution with that 
obtained from moments measured by ABN-2. 
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